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w This  study  was  conducted  to  determine  the  effect  that  the  proposed 
Marysville  Lake  pump-storage  hydropower  (2250-mw  power  plant)  project  would 
have  on  downstream  water  temperatures.  A one-dimensional  numerical  model  was 
used  for  simulation  and  prediction  of  temperatures  within  and  downstream  of 
Marysville  Lake.  Two  physical  hydraulic  models  were  used  for  study  and  de- 
scription of  the  hydrodynamic  response  of  the  project.  Information  from  the 
physical  models  was  necessary  to  determine  coefficients  used  by  algorithms 
within  the  mathematical  model.  The  mathematical  model  allowed  (Continued) 
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simulation  of  the  heat  exchange  characteristics  so  the  thermal  regimes 
within  and  downstream  of  the  lake  could  he  determined  for  various  hydrologic 
and  meteorologic  conditions  and  various  pumped-storage  hydropower  operations. 
Results  of  the  study  indicate  that  the  temperatures  should  he  within  the  ob- 
jective hand  desired  downstream  during  years  with  average  or  wetter  than 
average  hydrologic  conditions.  The  study  Indicated  that  with  the  ultimate 
2250-mw  power  plant  fall  temperattire  objectives  would  he  exceeded  by  a 
maximum  of  3°C  for  a U5-day  period  diiring  much  drier  than  average  years. 

With  the  initial  installed  capacity  of  1350  mw,  fall  temperatxire  objectives 
could  he  met  under  all  conditions  studied.. 
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PREFACE 


The  study  reported  herein  was  a continuation  of  a project  study 
that  was  authorized  by  the  Office,  Chief  of  Engineers  (OCE),  U.  S.  Army, 
on  l6  July  1975,  at  the  request  of  the  U.  S.  Army  Engineer  District, 

Sac  r ament  o ( SPK ) . 

This  study  was  conducted  to  determine  the  ability  of  the  proposed 
Marysville  Lake  Project  to  satisfy  downstream  temperature  objectives. 
Peaking  hydropower  and  pumped  storage  is  planned  for  Mairysville  Lake. 

The  power  producing  capacity  of  the  project  for  this  study  was  2250  mw. 

A previous  temperatiire  study  of  Marysville  Lake  was  conducted  by  the 
U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES)  and  reported  in 
WES  TR-H-7T-5,  Report  1.  That  study  was  concerned  with  the  ability  of 
the  Marysville  Lake  Project  designed  for  a 900-mw  capacity  to  satisfy 
downstream  temperature  objectives.  The  large  increase  in  flow  regimes 
brought  about  by  the  formulation  of  the  2250-mw  project  required  a re- 
assessment of  Marysville  Lake  thermal  characteristics  as  reported  herein. 

This  investigation  was  conducted  during  the  period  January  1977  to 
May  1977  in  the  Hydraxilics  Laboratory  of  WES  under  the  direction  of 
Messrs.  H.  B.  Simmons,  Chief  of  the  Hydraulics  Laboratory,  and  J.  L. 
Grace,  Jr.,  Chief  of  the  Structures  Division  and  Reservoir  Water 
Quality  Branch  (Physical).  The  study  was  conducted  by  Mr.  M.  S.  Dortch 
with  assistance  from  Messrs.  D.  G.  Fontane,  C.  H.  Tate,  Jr.,  and  D.  H. 
Merritt.  This  report  was  prepared  by  Mr.  Dortch  and  reviewed  by 
Mr.  Grace. 

Director  of  WES  during  this  study  was  COL  John  L.  Cannon,  CE. 
Technical  Director  was  Mr.  F.  R.  Brown. 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (Sl) 
UNITS  OF  MEASUREMENT 


U.  S.  customeu:y  units  of  measurement  used  in  this  report  can  be  con- 
verted to  metric  (Sl)  units  as  follows: 


Multiply 


feet 

miles  (U.  S.  statute) 
square  feet 
acre-feet 

cubic  feet  per  second 

feet  per  second  per 
second 


JJL 


0.301+8 

1.60931+1+ 

0.0929030I+ 

1233.1+82 

0.02831685 


To  Obtain 


O.30I+8 


metres 

kilometres 

square  metres 

cubic  metres 

cubic  metres  per 
second 

metres  per  second 
per  second 


Figure  1.  Location 
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Hydraulic  and  Mathematical  Model  Investigation 


PART  I : INTRODUCTION 

Pro.lect  Description 

1.  During  the  planning  stages  of  Marysville  Lake  (Figure  l),  the 
project  was  reformulated  such  that  the  hydropower  capacity  was  increased 
from  900  mw  to  an  ultimate  installation  of  2250  mw.  Most  of  the  project 
features  remained  unchanged  from  the  description  presented  in  the  pre- 
vious WES  report  (Report  l)*  that  pertained  to  the  900-mw  project.  The 
Marysville  Lake  pool  elevations  and  storage  capacity  did  not  change. 

The  increase  in  power  capacity  required  an  increase  in  discharge  rates 
between  the  forebay  and  afterbay  during  generation  and  pumpback.  Maxi- 
mum generation  flows  of  50,000  cfs**  for  the  900-mw  project  increased  to 
105,000  cfs  for  the  2250-mw  project.  Ptunpback  flow  rates  also  increased 
from  a maximum  of  9,000  cfs  (900-mw  project)  to  50,000  cfs  (2250-mw 
project).  Durations  of  flow  did  not  change  significantly;  thus  a greater 
volume  of  water  was  exchanged  during  generation  and  pumpba'rk  with  the 
2250-mw  project. 

2.  The  gross  storage  of  the  afterbay  was  increased  from 

lt0,U00  acre-feet  at  el  233t  to  80,U00  acre-feet  at  el  270.  However,  for 
the  study  years  presented  in  this  report,  a total  storage  of  only 
UU,300  acre-feet  was  required.  The  minimum  storage  of  the  afterbay  was 
increased  from  8,900  to  9,300  acre-feet. 

* D.  G.  Fontane  et  al.,  "Marysville  Lake  Hydrothermal  Study;  900-mw 
Project;  Hydraulic  and  Mathematical  Model  Investigation,"  Technical 
Report  H-7T-5,  Report  1,  May  1977,  U.  S.  Army  Engineer  Waterways 
Experiment  Station,  CE,  Vicksburg,  Miss. 

**  A table  of  factors  for  converting  U.  S.  customary  units  of  measure- 
ment to  metric  (SI)  units  is  presented  on  page  3. 
t All  elevations  (el)  cited  herein  are  in  feet  referred  to  mean  sea 
level . 
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3.  The  2250-mw  power  plant  would  consist  of  foiir  U50-mw  pump  tur- 
bines and  one  U50-mw  conventional  turbine.  Each  turbine  would  have  a 
selective-withdrawal  intake  system  consisting  of  vertically  spaced  ports 
opening  into  a common  wet  well.  Flow  would  pass  through  a single  se- 
lected port  (no  blending)  into  the  wet  well  and  from  the  wet  well  into 
the  penstock  intake.  Each  selective-withdrawal  port  would  have  an 
opening  25  ft  high  by  100  ft  wide.  Port  elevations  are  discussed  later 
in  this  report.  Based  on  recommendations  made  in  the  previous  study, 
it  was  decided  that  pumpback  flow  would  be  allowed  to  pass  only  through 
the  port  that  was  just  previously  used  for  withdrawal. 


Need  for  and  Purpose  of  Study 

4.  This  study  was  conducted  to  determine  the  effect  on  downstream 
water  temperatures  of  operating  the  Marysville  Lake  Project  with  a 
2250-mw  power  plant.  Temperature  regimes  of  peaking  hydropower,  pumped- 
storage  projects  can  be  highly  dependent  upon  pmpback  and  generation 
flows.  Because  of  the  much  larger  generation  and  piampback  flows  asso- 
ciated with  the  2250-mw  project,  it  was  necessary  to  completely  re- 


evaluate the  project. 


Approach 


5.  The  thermal  characteristics  within  and  downstream  of  Marysville 
Lake  were  evaluated  with  the  mathematical  simulation  model  that  was  used 
for  the  900-mw  study  described  in  Report  1.  With  greater  magnitudes  of 
flow  and  changes  in  the  operation  scheme  associated  with  the  2250-mw 
project,  it  was  anticipated  that  revisions  would  have  to  be  made  to  the 
mathematical  model  to  better  describe  the  hydrothermal  characteristics 

of  the  2250-mw  project. 

6.  As  in  the  900-mw  study,  physical  hydraulic  models  were  used  to 
obtain  an  improved  understanding  of  the  hydrodynamic  characteristics. 
Modifications  made  to  the  mathematical  model  were  based  on  physical 
model  resiilts.  Two  physical  models  were  used  in  the  study,  one  of  which 
was  also  used  during  the  900-mw  study. 
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PART  II:  PHYSICAL  MODEL  STUDIES 


Description 

Lake  model 

T.  The  distorted  scale  (l:l600  horizontal,  l:l60  vertical)  model 
of  Marysville  Lake  and  afterbay  that  was  used  for  the  900-mw  project 
study  was  also  used  for  the  2250-mw  study.  A complete  description  of 
the  model  is  presented  in  Report  1.  Because  the  geometry,  surface  area, 
and  storage-elevation  relation  of  Marysville  Lake  and  afterbay  did  not 
change  from  that  of  the  900-mw  project,  major  revision  of  the  model  was 
not  required.  However,  an  afterbay  storage  of  4i|,300  acre-feet  was 
needed  for  the  selected  study  years.  This  was  approximately  the  maximum 
storage  as  reproduced  in  the  existing  model  afterbay,  so  it  was  neces- 
sary to  add  some  additional  wall  height  to  the  afterbay  for  freeboard. 
Additionally,  the  intake-outlet  port  in  the  model  dam  had  to  be  enlarged 
to  account  for  the  larger  ports  in  the  intake  structiire  of  the  2250-mw 
project.  The  only  other  change  to  the  model  consisted  of  replacing  the 
reversible,  variable-speed  pump  with  a similar  pump  of  greater  capacity 
so  that  the  larger  generation  and  pumpback  flows  of  the  2250-mw  project 
could  be  simulated. 

Model  of  pimpback  jet 

8.  With  the  greater  volumes  of  water  pumped  back  each  day,  pump- 
back  mixing  associated  with  the  2250-mw  project  wo^lld  have  a greater 
effect  on  thermal  stratification  than  that  of  the  900-mw  project.  The 
algorithms  used  by  the  mathematical  model  to  describe  pumpback  mixing 
appropriately  have  a more  significant  influence  for  large  ptimpback  flows. 
Particularly,  the  gross  entrainment  coefficient,  E (reference  Report  l), 
has  a greater  influence  on  mixing  within  Marysville  Lake  for  the  2250-mw 
project  than  for  the  900-mw  project.  For  example,  when  the  p\mipback 
flow  is  9,000  cfs  and  E = 0.7  , the  entrained  flow  is  6,300  cfs.  For  a 
larger  piimped  flow  such  as  33,000  cfs  (typical  siunmer  condition, 

2250-raw  project),  an  entrained  flow  of  23,100  cfs  would  result  with 
E = O.T  . The  significance  of  the  value  of  a given  E is  much  greater 
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in  the  latter  case,  especially  when  considering  that  the  total  pumping 
time  and  the  lake  size  are  the  same  in  both  cases. 

9.  An  undistorted  1: 200-scale  model  of  a single  100-ft-wide  by 
25-ft-high  port  was  used  to  study  pumped  flow  and  evaluate  E . The 
port  was  located  within  the  end  wall  of  an  existing  3-ft-wide  by  2-ft- 
deep  by  25-ft-long  rectangular  transparent  plastic  flume.  With  the 
bottom  of  the  flume  corresponding  to  prototype  el  200  (approximate  res- 

f ervoir  bottom,  Marysville  Lake),  the  port  was  fixed  at  1.50  ft  above 

the  bottom  or  prototype  el  500.  Water  was  pumped  from  a supply  tank, 
through  a rotameter  and  valve,  through  the  wet  well  and  port,  into  the 
; flume.  With  the  inside  of  the  flxmie  end  wall  representing  the  upstream 

'i  face  of  the  dam,  the  wet  well  was  attached  on  the  outside  of  the  end 

wall  covering  the  port.  The  model  wet  well  reproduced  the  preliminary 
wet  well  design  for  the  prototype,  thus  providing  a reasonable  velocity 
distribution  through  the  port.  The  pumped  flow  was  measured  with  the 
rotameter  and  scaled  to  prototype  flow  by  means  of  Froude  scaling 
criteria. 

Tests  and  Results 

Entrainment  characteristics 

10.  Steady-state,  short-term  pumpback  flow  was  studied  with  the 
\indistorted  model  to  determine  the  entrainment  characteristics.  Density 
stratification  resulting  from  thermal  stratification  within  Marysville 
Lake  was  simulated  in  the  fl\jme  with  saline  and  fresh  waters.  The 
density  of  the  pumpback  water  was  varied  from  test  to  test,  but  was 
always  greater  than  the  density  of  the  surface  water  and  less  than  the 
density  of  bottom  water  of  the  pool.  A constant  pumpback  flow  of 
l6,800  cfs  was  simulated  in  the  model.  This  flow  is  typical  of  the 
pumped  flow  passing  through  a single  port  from  a single  pump  turbine. 

An  average  summer  pool  of  el  535  was  used  throughout  the  tests. 

11.  A density  probe  consisting  of  a conductivity  and  temperature 
sensor  calibrated  against  known  specific  gravity  was  used  to  obtain 
density  measurements.  Flow  measurements  in  the  pool  were  obtained  from 
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dye  streak  displacements.  Dye  particles  were  dropped  into  the  pool.  As 
the  particles  fell,  vertical  dye  streaks  were  created.  Fluid  motion 
created  a dye  streak  displacement.  By  recording  the  dye  streak  dis- 
placement with  video  equipment,  velocities  could  be  obtained  during 
playback.  By  integrating  the  velocity  distributions  with  respect  to 
depth,  unit  discharges  were  calculated.  Practically  uniform  flow  was 
observed  across  the  width  of  the  flume  at  stations  outside  of  the  spread- 
ing and  turbulent  mixing  zone  of  the  Jet;  therefore  the  unit  discharges 
were  multiplied  by  the  width  to  obtain  volume  flow  rates.  The  velocity 
profiles  were  taken  outside  of  the  turbulent  spreading  and  mixing  zone 
of  the  pumpback  jet.  It  was  observed  that  the  lateral  spread  of  the  jet 
was  not  restricted  by  the  flume  walls. 

12.  With  a knowledge  of  the  volume  flow  rates  of  the  pumpback 
current  and  the  pumping  rate,  it  was  possible  to  estimate  the  total 
amount  of  flow  that  was  entrained  by  and  mixed  with  the  pmpback  flow. 

The  entrainment  of  flow  causes  the  flow  of  the  pumpback  current  to  be 
greater  than  the  flow  rate  of  the  pumpback  jet.  The  pumpback  jet  actu- 
ally induces  flow.  The  volume  rate  of  flow  of  the  pumpback  current  can 
be  described  by 


Q = Q + EQ 
CO  o 


where 

= volume  flow  rate  of  the  ptunpback  current,  cfs 
= volijme  flow  rate  of  the  pumpback  jet,  cfs 
E = total  entrainment  coefficient 

It  is  necessary  to  know  the  entrainment  coefficient  in  order  to  properly 
acco\int  for  the  budgeting  of  water  and  temperature  within  the  pool. 

13.  Most  of  the  entrained  flow  was  observed  to  come  from  the 
layers  of  the  pool  corresponding  to  the  elevation  of  the  port  through 
which  the  jet  entered  the  pool.  This  information  is  needed  for  mathe- 
matical model  input  and  is  mentioned  again  in  the  section  on  mathemati- 
cal model  modifications. 

lU.  There  was  some  difficulty  in  obtaining  values  of  the 
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entrainment  coefficient,  E . Because  the  pumpback  flow  causes  distur- 
bances (water  displacement)  in  the  pool,  internal  gravity  waves  are 
created  that  complicate  measurement  of  the  advected  flow.  The  internal 
waves  caused  displacements  of  the  dye  streaks,  making  it  diffic\J.t  to 
distinguish  whether  the  motion  was  due  to  advection  or  waves.  Numerous 
dye  streak  displacements  were  recorded  during  each  test  so  that  several 
velocity  profiles  could  be  obtained  and  averaged  in  an  effort  to  damp 
periodic  motion  due  to  the  internal  waves.  The  average  velocity  profile 
was  used  to  obtain  a value  for  E for  each  test.  An  average  value  for 
E of  about  1.0  was  obtained  from  the  tests.  For  mathematical  model 
purposes,  a value  of  E = 1.0  was  applied  to  the  flow  from  the  pump 
turbine. 

15.  To  substantiate  the  above  test  resiilt,  another  procedure  was 
used  to  estimate  E . As  presented  in  Report  1,  there  is  a relation 
that  fairly  well  describes  the  thickness  of  an  interflow,  such  as  that 
due  to  a pumpback  current.  The  thickness,  D , can  be  calculated  from 

D = It.l 

where 

W = average  reservoir  width  at  the  elevation  of  the  pimipback 
current , ft 

Ap  = density  difference  of  the  epilimnion  and  hypolimnion,  g/cc 

p = average  density  of  the  ptimpback  current,  g/cc 

''  2 

g = acceleration  due  to  gravity,  ft/sec 

By  dying  the  pimped  flow,  it  was  easy  to  determine  the  thickness,  D , 
of  the  pumpback  current.  Knowing  Ap  and  W , could  be  calculated 

from  Equation  2 and  substituted  into  Equation  1 to  solve  for  E . Values 
of  E obtained  this  way  compared  favorably  with  those  obtained  from  the 
average  velocity  computations. 

Mixing  due  to  pumpback 

16.  As  discussed  in  Report  1,  the  pumpback  current  alters  the 
density  structure  as  it  passes  through  the  interflow  zone.  To  account 
for  the  displacement  and  mixing  caused  by  the  pumpback  current,  a mixing 
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technique  that  was  used  in  the  mathematical  model  for  the  900-mw  study 
was  also  used  for  this  study.  The  technique,  which  is  discussed  in 
greater  detail  in  Report  1,  is  based  on  the  concept  of  a portion  of  each 
of  the  layers  within  the  zone  of  the  pumphack  current  being  removed, 
mixed  together,  returning  the  mixed  water  to  the  layers  within  the  zone, 
and  then  mixing  within  each  layer.  The  portion  of  water  removed  and 
mixed  from  each  layer  is  computed  by  multiplying  a mixing  coefficient 
times  the  volume  of  the  layer.  The  mixing  technique  is  applied  once 
during  each  day  of  mathematical  model  simulation.  Physical  model  data 
indicated  that  the  mixing  coefficients  decay  exponentially  with  vertical 
distance  above  or  below  the  ptimpback  ciirrent  inflow  layer.  The  pumphack 
current  inflow  layer  is  the  layer  of  the  pool  where  the  lake  water 
density  most  nearly  equals  the  average  density  of  the  pumphack  ciirrent, 

. The  form  of  this  exponential  equation  is 

-BX. 

= Ae  ^ (3) 

where 

= mixing  coefficient 
e = natural  logarithmic  base  (2.7183) 

A,B  = coefficients 

= distance  from  the  pumphack  current  inflow  layer  to  layer 
i = layer  number 

IT.  Tests  were  conducted  with  the  distorted-scale  model  of 
Marysville  Lake  to  evaluate  the  coefficients  A and  B for  Equation  3. 
Generation  and  pumphack  operations  were  simulated  for  conditions  cor- 
responding to  an  average  hydrologic  year.  Density  profiles  were  ob- 
tained from  the  model  periodically.  Through  analysis  of  changes  in  the 
density  profiles  over  known  time  periods,  it  was  possible  to  evaluate 
. A fit  of  the  data  provided  the  coefficients  A and  B used  in 
Equation  3.  The  values  of  A and  B that  were  found  to  be  appropriate  for 
the  2250-mw  project  were  0.l6U  and  0.035,  respectively. 

Afterbay  mixing 

18.  Although  the  afterbay  can  be  assumed  to  be  fully  mixed  for 


mathematical  modeling  (Report  l),  it  is  necessary  to  properly  budget  the 
advection  of  heat  to  and  from  the  afterbay.  As  discussed  in  Report  1, 
this  is  accomplished  by  accounting  for  the  pregeneration  afterbay  water 
(fully  mixed)  and  the  Marysville  Lake  generation  water  that  is  pumped 
back  in  a day.  Pregeneration  afterbay  water  and  Marysville  Lake  genera- 
tion water  that  is  not  pumped  back  is  mixed  in  the  afterbay  to  obtain  a 
volume-weighted  average  temperature  of  afterbay  water  after  generation. 

19.  Tests  were  conducted  with  the  physical  lake  model  to  de- 
termine how  much  pregeneration  afterbay  water  constituted  the  pumpback. 
With  this  knowledge,  mixing  of  pregeneration  afterbay  water  with  gener- 
ation water  could  be  properly  accotinted  for.  Additionally,  this  infor- 
mation allows  the  average  temperattire  of  the  pumpback  water  to  be  cal- 
culated. As  in  the  900-mw  study,  fluorescent  dye  was  used  to  trace  the 
pregeneration  afterbay  water  being  pumped  back  for  various  operation 
conditions.  These  testing  procedures  are  explained  in  Report  1. 

20.  The  ratio  of  V , the  total  volme  of  water  pumped  back  in 

a day,  to  , the  total  volume  of  water  released  diiring  generation  in 
a day,  was  found  to  be  a predominant  factor  influencing  the  amount  of 
pregeneration  afterbay  water  comprising  the  pumpback.  Additionally, 
the  ratio  of  , the  pregeneration  afterbay  volume,  to  V was  found 
to  be  an  important  parameter  for  determining  the  pumpback  constituents. 
From  the  fluorescent  dye  tests,  , the  volume  of  pregeneration  after- 
bay water  that  is  pumped  back  in  a day,  was  determined  for  various  op- 
eration conditions.  Results  of  these  tests  are  shown  by  a plot  of  the 
parameters  versus  ^ (Plate  1).  A least  squares 

fit  of  these  two  parameters  resulted  in  the  equation 

/v  v_ 

Test  conditions  shown  in  Plate  1 were  typical  of  most  of  the  conditions 
encountered  during  the  mathematical  simulations.  Because  the  parameter 
V^/Vp  is  undefined  by  test  results  for  values  of 

than  0.50,  V./V  was  assumed  constant  at  O.UO  in  the  mathematical  model 
A p 

for  values  of  ^^p/Vg^  greater  than  0.50. 
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PART  III:  MATHEMATICAL  MODEL 


I 


Model  Modifications 


21.  Modifications  made  to  the  mathematical  model  used  in  the 
900-mw  study  consisted  of  incorporating  the  physical  model  results  that 
were  discussed  in  the  previous  section.  A complete  description  of  the 
mathematical  model  can  be  found  in  Report  1.  The  modifications  con- 
sisted of  the  following: 

a.  Changed  the  value  of  the  entrainment  coefficient,  E , 
from  O.T  (900-mw  study)  to  1.0. 

b.  Changed  the  coefficients  of  the  pumpback  mixing  equation 
(Equation  3)  from  A = 0.16T5  and  B = 0.1l8  (900-mw 
study)  to  A = 0.l6U  and  B = 0.035. 

c_.  Incorporated  a new  equation  (Equation  it)  to  describe 
mixing  in  the  afterbay. 

d.  Modified  the  description  of  entrained  flow  resulting  from 
pumpback  to  allow  entrainment  from  several  layers  of  the 
pool. 

The  entrainment  layers  consist  of  all  layers  between  the  elevations  of 
the  invert  and  top  of  the  port.  Since  the  ports  are  25  ft  high  and  the 
mathematical  layer  size  is  5 ft,  there  are  five  pmpback  entrainment 
layers.  The  five  layers  contribute  equally  to  the  total  entrained  flow. 


22.  The  general  data  changes  Liaue  for  this  study  consisted  of: 
The  number  of  power  units  was  set  equal  to  5. 

2 

b.  The  individual  port  area  was  changed  to  2500  ft  . 

c_.  The  minimum  afterbay  volume  was  fixed  at  9300  acre-feet. 

An  isothermal  starting  temperature  of  12°C  was  used  for 
the  simulations  rather  than  8°C,  which  was  used  in  the 
900-raw  study. 

e^.  The  hydrology  was  revised  to  reflect  partial  impairment 
usage  of  upstream  water  rights. 

A starting  temperature  of  12®C  was  fo\ind  to  be  more  appropriate  for 
simulating  this  project.  The  use  of  an  8°C  starting  temperature  was 
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based  on  observed  temperature  profiles  taken  at  nearby  Englebright  and 
New  Bullards  Bar  Reservoirs,  which  are  practically  isothermal  at  about 
8°C  at  the  beginning  and  end  of  a calendar  year.  Both  of  these  lakes 
have  high  level  releases.  When  Marysville  Lake  is  simulated  without 
pumped  storage  and  with  a high-level  release  (all  year),  a starting  tem- 
perature of  8°C  should  be  used  because  the  lake  temperatures  return  to 
about  8°C  or  9°C  isothermal  at  the  end  of  the  year.  However,  with  pumped 
storage  and  low-level  releases  in  the  fall,  the  lake  temperatures  return 
to  only  about  12°C.  This  is  reasonable  since  mixing  resulting  from 
pumped-storage  causes  more  heat  to  enter  the  pool,  and  low-level  releases 
in  the  fall  do  not  allow  the  advection  of  heat  out  of  the  lake  from  the 
warmer  upper  layers.  Simulations  were  conducted  with  a starting  temper- 
a’:ure  of  8°C  and  12°C.  Temperatures  predicted  by  the  two  simulations 
differed  during  the  winter  and  early  spring  months;  but  after  the  model 
was  allowed  to  simulate  several  spring  months,  almost  identical  tempera- 
ture regimes  were  found. 

23.  The  number  of  selective-withdrawal  intakes  for  each  power 
\mit  was  reduced  from  6 for  the  study  to  5 for  the  2250-mw  study. 

A conclusion  of  the  previous  study  was  that  the  number  of  withdrawal 
elevations  could  be  reduced  without  adversely  affecting  the  desired  re- 
lease temperatures.  The  elevations  of  several  ports  also  were  changed 
from  the  previous  study.  The  port  locations  used  for  this  study  are 
shown  below: 

Center-line  Elevation 
Port  ft  msl 

1 530.0 

2 1*90.0 

3 1*1+0. 0 

1*  370.0 

5 252.5 


The  port  size  used  for  this  study  was  100  ft  wide  and  25  ft  high. 
2l*.  The  study  years  used  for  this  study  consist  of: 


ll* 


m 


1 

J 


Study  year  1 

1962 

Hydrology  (average) 

1962 

Meteorology  (average) 

Study  year  2 

1931 

Hydrology  (dry) 

1967 

Meteorology  (hot) 

Study  year  3 

191*2 

Hydrology  (wet) 

1963 

Meteorology  (cold) 

study  years  1 and  3 are  the  same  hydrologic  and  meteorologic  years  as 
those  used  for  the  900-mw  study.  For  study  year  2,  hydrologic  year  1931 
(a  dry  year  with  pumphack)  was  substituted  for  hydrologic  year  193^ 

(a  dry  year  without  piompback).  The  meteorologic  year  1967  was  also  used 
for  study  year  2 in  the  900-mw  study.  The  hydrologic  routings  supplied 
by  SPK  for  the  900-mw  study  were  obtained  for  the  condition  of  iMll  im- 
pairments of  water  rights.  Full  impairments  mean  that  the  preimpoond- 
ment  water  rights  upstream  of  Marysville  Lake  are  assumed  to  be  fully 
utilized  throughout  each  year.  However,  the  water  rights  are  not  fully 
utilized  (partial  impairments)  during  the  winter  months  and  portions  of 
the  spring  and  fall.  This  results  in  largei'  inflows  and  outflows  than 
would  occur  with  full  impairments.  Partial  impairment  of  water  utili- 
zation is  a more  realistic  description  of  the  hydrology  and  was  used  by 
SPK  to  develop  the  hydrologic  routings  and  project  operations  for  the 
2250-mw  study.  Inflows  to  Marysville  Lake,  afterbay  releases,  pool 
storage,  and  irrigation  diversions  at  Daguerre  Point  for  study  years  1, 
2,  and  3 are  presented  in  Tables  1-3.  The  daily  operation  schedule 
planned  for  the  2250-mw  project  is  quite  different,  in  terms  of  flow 
rates,  from  that  planned  for  the  900-mw  project.  The  2250-mw  operation 
schedules  provided  by  SPK  for  the  three  selected  study  years  are  pre- 
sented in  Tables  1+-6.  The  flow  rates  and  durations  shown  in  Tables  U-6 
are  monthly  average  daily  values. 

Model  Output 

25.  Mathematical  simulations  of  the  2250-mw  project  were  con- 
ducted for  all  three  study  years.  Predicted  in-lake  temperature  pro- 
files for  the  three  study  years  are  presented  in  Plates  2,  U,  and  6. 
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Predicted  release,  afterbay,  and  downstream  temperatiires  for  the  three 
study  years  are- presented  in  Plates  3,  5,  and  7.  The  selective- 
withdrawal  port  selection  throughout  each  year  is  also  indicated  in 
Plates  3,  5,  and  7.  Temperature  objectives  are  to  be  satisfied  "down- 
stream” at  the  confluence  of  the  Yuba  and  Feather  Rivers.  These  ob- 
jectives are  the  same  as  those  used  for  the  900-mw  study  and  are  defined 
by  the  downstreeun  objective  band  shown  in  Plates  3,  5,  and  7. 
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PART  IV:  DISCUSSION  OF  RESULTS 

26.  Downstream  temperat\ires  stay  within  the  objective  temperat\ire 
band  for  study  years  1 and  3 as  shown  by  Plates  3 and  T.  Through  care- 
f\al  port  selection,  it  should  be  possible  to  stay  within  the  band  for 
years  similar  to  study  years  1 and  3.  For  prototype  operations,  it 
would  be  advisable  to  select  ports  based  on  the  results  of  a mathe- 
matical operations  model  of  the  project. 

27.  It  was  not  possible  to  operate  the  2250-mw  project  as  pres- 
ently planned  and  stay  within  the  desired  temperature  band  at  all  times 
for  study  year  2 (Plate  5)«  Should  the  2250-mw  project  be  operated  with 
ptuupback  during  hot,  dry  (low  pool)  years,  difficulty  in  staying  within 
the  band  could  occur.  The  release  temperatxires  of  study  year  2 can  be 
held  within  the  band  diiring  the  summer,  but  this  will  result  in  even 
warmer  releases  in  the  fall.  The  warmer-than-desired  summer  tempera- 
tures (Plate  5)  are  caused  by  the  hotter-than-average  meteorologic  con- 
ditions that  create  additional  warming  within  the  stream  between  the 
afterbay  and  the  confluence  of  the  Yuba  and  Feather  Rivers.  The  warmer- 
than-desired  fall  releases  are  due  to  the  shortage  of  cold  water  in  the 
Marysville  pool.  The  combination  of  a dry  year  and  pumpback  reduced 
the  availability  of  cold  water.  A dry  year  without  pumpback  co\ild  prob- 
ably meet  fall  temperature  objectives  as  was  demonstrated  by  study 

year  2 of  the  900-mw  study.  A hot  year  with  pumpback  and  average  hy- 
drologic conditions  can  also  meet  fall  temperature  objectives  as  shown 
by  Plate  8.  The  release  plots  in  Plate  9 show  that  for  a dry  year  with 
average  meteorologic  conditions,  the  fall  objectives  still  could  not  be 
met.  The  downstream  temperatures  were  cooler  during  the  svunmer  than 
those  shown  in  Plate  5;  the  objective  band  could  be  met  during  the 
summer,  but  warmer  fall  releases  woixld  result.  Therefore,  the  critical 
factor  for  study  year  2 is  that  it  is  dry.  Because  it  is  dry  the 
Marysville  pool  is  lower  than  usual,  thus  containing  less  cold-water 
storage.  Mixing  that  results  from  pumpback  causes  warming  of  the  cold 
water  that  is  available. 

28.  A simulation  was  conducted  to  determine  how  the  downstream 
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water  temperatures  predicted  for  study  year  2 with  the  project  differed 
from  those  predicted  for  the  same  study  year  without  the  project.  To 
make  this  comparison,  predicted  downstream  temperatures  for  study  year  2 
were  plotted  with  predicted  natiiral  downstream  temperatures  (without  the 
project)  for  study  year  2.  The  natxiral  downstream  temperatiares  were 
developed  hy  routing  inflows  to  Marysville  from  the  headwater  region  of 
Marysville  Lake  to  the  confluence  of  the  Yuba  and  Feather  Rivers.  Heat 
exchange  for  this  stream  routing  was  applied  in  the  same  manner  as  that 
used  to  obtain  the  downstream  temperatures  for  the  lake  simulations  and 
is  explained  in  Report  1.  As  shown  by  Plate  10,  the  predicted  down- 
stream temperat\ires  with  and  without  the  project  are  very  similar  ex- 
cept in  the  fall  where  the  temperatures  with  the  project  are  warmer. 

29.  The  results  discussed  thus  far  are  associated  with  the 
2250-mw  project.  As  presently  planned,  the  project  would  have  an 
initial  installed  capacity  of  1350  mw,  consisting  of  one  ii50-mw  con- 
ventional turbine  and  two  U50-mw  pump  turbines.  Design  of  the  power 
plant  would  provide  for  future  installation  of  two  additional  U50-mw 
pump  turbines  when  the  future  need  for  power  is  confirmed,  bringing 
the  total  capacity  of  the  plant  to  2250  mw.  Because  there  is  less 
pumpback,  it  would  be  less  difficult  to  meet  downstream  temperature 
objectives  with  the  1350-mw  power  plant  than  with  the  2250-mw  power 
plant.  Because  temperature  objectives  could  be  met  for  study  years  1 
and  3 with  the  2250-mw  power  plaint,  it  should  be  even  less  difficult 
to  meet  objectives  during  these  study  years  with  the  1350-mw  power 
plant.  However,  because  the  downstream  temperatures  were  outside  the 
objective  beuid  during  part  of  study  year  2 (2250  mw),  it  was  not  known 
whether  the  temperatures  could  be  maintained  within  the  band  during 
this  study  year  with  the  1350-mw  project. 

30.  Simulations  therefore  were  conducted  to  determine  the  effect  j 

that  the  project  with  a 1350-mw  power  plant  would  have  on  downstream  | 

water  temperatures  during  study  year  2 (dry,  hot  year).  Coefficients 

used  by  the  mathematical  model  for  the  2250-mw  simulations  were  also 
used  for  these  simulations.  This  is  considered  to  be  a conservative 
assumption.  The  hydrologic  input  is  shown  in  Table  2.  Operation 
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schedules  are  presented  in  Table  7-  Restilts  of  these  simiilations  are 
presented  in  Plates  11  and  12.  As  shown  by  Plate  12,  downstream  tem- 
peratures are  within  the  objective  band  dixring  the  fall.  Downstream 
temperatures  were  as  much  as  2°C  above  the  band  from  late  July  to 
mid-September  because  of  warming  that  occ\ars  between  the  afterbay  and 
the  confluence  of  the  Yuba  and  Feather  Rivers.  By  releasing  colder 
water  in  siimmer,  objective  temperatures  could  be  maintained  in  the 
summer,  but  this  wovild  deplete  cold-water  storage  to  the  extent  that 
warmer-than-de sired  releases  would  restilt  in  the  fall. 


PART  V:  SUMMARY 


31.  This  study  indicated  that  water  temperatures  downstream  from 
Marysville  Lake  project  could  be  maintained  within  the  desired  tempera- 
ture range  during  years  with  average  or  wetter-than-average  hydrologic 
conditions,  with  either  a 1350-mw  or  a 2250-mw  pumped-storage  puv^er 
plant  operating.  The  model  studies  indicated,  however,  that  with  the 
2250-mw  project  operating  under  the  much-drier-than-average  conditions 
experienced  in  1931,  the  desired  downstream  water  temperatures  would  he 
exceeded  by  a maximum  of  3°C  from  mid-July  through  mid-September  and  by 
about  2°C  in  mid-October,  gradually  decreasing  to  desired  levels  by  the 
end  of  November.  With  the  initially  installed  1350-mw  power  plant  oper- 
ating under  1931  conditions,  desired  temperatures  would  be  met  during 
the  fall,  but  would  be  exceeded  by  as  much  as  2°C  from  late  July  through 
mid-September.  By  modifying  the  amount  of  pumpback  during  critical  pe- 
riods or  by  taking  other  measures,  it  is  probable  that  the  desired  water 
temperatures  could  be  maintained  throughout  the  year  under  1931  condi- 
tions. Other  measures  that  could  be  taken,  for  instance,  might  involve 
release  of  cold  water  that  would  be  available  within  the  Dry  Creek  Arm 
(see  Report  1,  paragraph  33)  down  Dry  Creek  to  mix  with  the  water  re- 
leased from  the  afterbay  on  the  Yuba  River  or  to  excavate  the  connecting 
channel  deeper  to  allow  the  cold  water  below  the  connecting  channel  in- 
vert in  the  Dry  Creek  Arm  to  mix  with  water  in  the  Yuba  River  Arm  (main 
part  of  the  lake). 
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I Table  1 

J Hydrologic  Input  Data.  Partial  Impairments 

Hydrologic  Year  1962 


Irrigation 
Withdrawal 

I End  of  Month  at 


Month 

Inflow 
cfs  per  day 

End  of  Month 
Lake  Storage 

X 1000  AF* 

Lake 

Elevation 

ft/msl 

Afterbay 
Release 
cfs  per  day 

Daguerre 

Point 

cfs  per  day 

Jan 

1563 

676.9 

521 

15I+O 

0 

Feb 

k331 

690.7 

523 

lt080 

0 

Mar 

2717 

773.1 

537 

1380 

50 

Apr 

2501 

772.0 

537 

2520 

310 

May 

22kk 

769.1 

537 

2290 

800 

Jun 

236h 

765.6 

536 

2U20 

1000 

Jul 

2288 

760.9 

535 

2360 

1130 

Aug 

177i+ 

756.8 

535 

l81t0 

1030 

Sep 

128U 

71+9.8 

53!+ 

lUOO 

800 

Oct 

51I+2 

736.0 

531 

5370 

360 

Nov 

12U5 

676.0 

521 

2250 

0 

Dec 

2223 

676.0 

521 

2220 

0 

Table  2 

Hydrologic  Input  Data,  Partial  Impairments 
Hydrologic  Year  1931 


Month 

Inflow 
cfs  per  day 

End  of  Month 
Lake  Storage 

X 1000  AF» 

End  of  Month 
Lake 

Elevation 

ft/msl 

Afterbay 
Release 
cfs  per  day 

Irrigation 

Withdrawal 

at 

Daguerre 

Point 

cfs  per  day 

Jan 

1669 

677.1 

521 

1660 

0 

Feb 

1830 

677.6 

521 

1820 

0 

Mar 

2017 

677:6 

521 

2020 

50 

Apr 

lh62 

676.5 

521 

llt80 

310 

May 

1306 

673.0 

520 

1360 

800 

Jun 

1281 

669.7 

520 

1330 

QUO 

Jul 

112i^ 

655.1 

517 

1350 

930 

Aug 

1039 

639.6 

511+ 

1280 

860 

Sep 

81T 

622.9 

511 

1090 

670 

Oct 

998 

621.  U 

511 

1020 

3h0 

Nov 

857 

621.  U 

511 

860 

0 

Dec 

1311 

667.5 

519 

560 

0 

AF  = acre-feet 


Table  3 

Hydrologic  Input  Data,  Partial  Impairments 
Hydrologic  Year  19^2 


Irrigation 

Withdrawal 

End  of  Month  at 


Month 

Inflow 
cfs  per  day 

End  of  Month 
Lake  Storage 

X 1000  AF* 

Lake 

Elevation 

ft/msl 

Afterbay 
Release 
cfs  per  day 

Daguerre 

Point 

cfs  per  day 

Jan 

6801 

686. 5 

523 

6630 

0 

Feb 

77lU 

706.6 

526 

7360 

0 

Mar 

2909 

795.3 

5^1 

IU70 

50 

Apr 

532it 

89i*.Jt 

557 

3570 

310 

May 

U956 

916.0 

560 

1*690 

800 

Jun 

3788 

916.0 

560 

3780 

1000 

Jul 

3091 

910.8 

559 

3170 

1130 

Aug 

1717 

906.2 

559 

1790 

1030 

Sep 

1258 

800.0 

5I+2 

30U0 

800 

Oct 

1258 

736.0 

531 

2290 

360 

Nov 

973 

676.0 

521 

1980 

0 

Dec 

233h 

676.0 

521 

2330 

0 

* AF  = acre-feet. 
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leration  Schedule.  Pai'tial  Impairments 
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Note:  No  generation  or  pcampback  on  Sundays.  There  are  two  pumpback  periods  during  a day  for  several  months. 
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Note:  No  generation  or  ptunpback  on  Sundays.  There  are  two  pumpback  periods  during  a day  for  several  months. 
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Note:  No  generation  or  pumpback  on  Sundays.  There  are  two  pumpback  periods  during  a day  for  severed,  months. 
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